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ABSTRACT 

The transcription factor Msn2 mediates a significant 
proportion of the environmental stress response, 
in which a common cohort of genes changes ex- 
pression in a stereotypic fashion upon exposure 
to any of a wide variety of stresses. We have ap- 
plied genome-wide chromatin immunoprecipitation 
and nucleosome profiling to determine where Msn2 
binds under stressful conditions and how that bind- 
ing affects, and is affected by, nucleosome posi- 
tioning. We concurrently determined the effect of 
Msn2 activity on gene expression following stress 
and demonstrated that Msn2 stimulates both activa- 
tion and repression. We found that some genes re- 
sponded to both intermittent and continuous Msn2 
nuclear occupancy while others responded only to 
continuous occupancy. Finally, these studies docu- 
ment a dynamic interplay between nucleosomes and 
Msn2 such that nucleosomes can restrict access of 
Msn2 to its canonical binding sites while Msn2 can 
promote reposition, expulsion and recruitment of nu- 
cleosomes to alter gene expression. This interplay 
may allow the cell to discriminate between different 
types of stress signaling. 

INTRODUCTION 

Regulation of eukaryotic gene expression involves a com- 
plex interplay among transcription factors, core transcrip- 
tional machinery and the chromatin template on which 
these factors operate. A number of studies over the last sev- 



eral years have documented that the chromatin structure 
across a cell's genome remains well defined and remarkably 
static under all conditions (1-3). Generally, well-positioned 
nucleosomes bracket the promoter region of most genes to 
maintain a nucleosome-depleted region (NDR) upstream of 
the transcriptional start site of the gene, with nucleosomes 
assuming a well-ordered periodic array extending into the 
coding region with periodicity diminishing with increasing 
distance from the promoter (1,4-6). This chromatin struc- 
ture serves an instructive role in transcription factor bind- 
ing, with factors able to bind to their cognate sites lying 
within the NDR but unable to bind to those sites occluded 
by nucleosomes in other regions (3,7,8). Against this back- 
drop of static chromatin structure, nucleosome depletion 
around the NDR is in some cases associated with transcrip- 
tional activation and nucleosome recruitment to the NDR 
associated with transcriptional repression (9-12). This lo- 
cal reorganization depends on the action of chromatin re- 
modeling factors that slide, evict or recruit nucleosomes 
(2,11,13-19). These rearrangements also occur in concert 
with transcription factor binding and transcriptional repro- 
gramming, although the causal nature of those relations is 
not entirely clear. To address this question, we have exam- 
ined transcriptional reprogramming and nucleosome rear- 
rangements associated with the yeast stress response. 

All cells mount a rapid adaptive response to a new and 
stressful environment and that response generally includes 
substantial transcriptional reprogramming. The transcrip- 
tional response of yeast cells to any of a wide variety of 
stresses, including heat shock, oxidative agents, nutrient de- 
pletion and hypo- and hyperosmolarity, comprises a stereo- 
typic repression and induction of the same large number of 
genes independent of the particular type of stress, referred 
to as the environmental stress response (ESR), as well as 
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changes specific to the individual stress (20,21). This tran- 
scriptional reprogramming has minimal impact on the sur- 
vival of cells to the eliciting stressful condition but does 
serve to protect cells to subsequent stresses (22-25). 

The structurally related, stress-responsive transcription 
factors, Msn2 and Msn4, mediate a major component of 
the ESR (21,26). These two transcription factors reside 
in the cytoplasm in unstressed cells, due to active export 
from the nucleus by the Msn5 exportin machinery and to 
restricted import as a consequence of protein kinase A 
(PKA) catalyzed phosphorylation and inhibition of the nu- 
clear import signals on the proteins (27-30). A number of 
microfluidics-based single-cell time lapse studies have doc- 
umented that acute stress causes rapid cycling of Msn2 
and Msn4 into and out of the nucleus, due to inhibition 
of PKA and activation of protein phosphatase 1 and the 
Snfl adenosine-monophosphate (AMP)-activated kinase 
(31-35). Cells exhibit idiosyncratic patterns of Msn2 nu- 
clear cycling such that genetically identical cells under the 
same stress condition show markedly different patterns of 
cycling. Moreover, different stresses elicit different classes 
of nuclear localization patterns (34,35). How Msn2 cycling 
relates to the transcriptional output from Msn2 remains to 
be resolved, although recent results suggest that different 
promoters respond to Msn2 cycling in different ways (36). 

Once in the nucleus, Msn2 can bind to stress response el- 
ements (STREs) within the genome to alter transcription 
of genes neighboring the sites (26,37). Approximately 3000 
STREs reside upstream of yeast genes, but only a fraction 
of these serve as binding sites for Msn2. Many of these sites 
are likely occluded by positioned nucleosomes that prevent 
access to Msn2 (38). Moreover, since each cell contains only 
100-200 Msn2 molecules (39), formation of stable Msn2 
complexes with a large number of STREs within a single 
cell's genome would not be possible. One explanation for 
the rapid dynamics of Msn2 localization may be to facilitate 
sampling of multiple sites by individual Msn2 molecules. 
Whether different stresses affect the selection of different 
subsets of sites — either by modifying Msn2's deoxyribonu- 
cleic acid (DNA) binding recognition region or by altering 
the accessibility of different sites — has not been extensively 
explored. 

Several previous studies have addressed the localization 
of Msn2 binding on a genome-wide basis. Venters et al. (40) 
examined global Msn2 binding in response to heat shock 
in the context of a much larger study to map the major- 
ity of transcription factor binding sites in yeast. More re- 
cently, Huebert et al. (41) mapped the location of Msn2 
binding over time over the entire genome following treat- 
ment of cell with peroxide and correlated that binding with 
genome-wide changes in nucleosome positioning. Here we 
examine the binding of Msn2 to genomic sites in response 
to a nutritional stress. We also measure the global nucleo- 
some architecture before and after application of the stress 
both in the presence and absence of Msn2 in order to ad- 
dress the extent to which Msn2 binding influences and is in- 
fluenced by nucleosomes. Finally, we assess the sufficiency 
and necessity of Msn2 binding on changes in expression of 
each associated gene to determine the effect on transcrip- 
tion elicited by Msn2 binding. Our results document an ex- 
tensive interplay between nucleosome binding and Msn2 



binding such that nucleosome occlusion can restrict Msn2 
binding in some conditions but in other cases Msn2 binding 
leads to repositioning of occluding nucleosomes. These re- 
sults point to a complex choreography between general and 
specific transcription factors in order to mount a coherent 
transcriptional program. In a companion paper (Elfving et 
al. 2014, submitted), we also examine the role of Mediator 
in this process. 



MATERIALS AND METHODS 

Yeast strain growth and construction 

Standard methods were used to grow, maintain and con- 
struct strains (42). Strains used in this study are listed in 
Table 1. 

The myc-tagged Msn2 strain was made by replacing the 
stop codon in the corresponding ORF (YMR037C) with a 
G8-13xMyc-KanMx6 construct. This construct was made 
by polymerase chain reaction (PCR) from the template plas- 
mid pFA6a-13Myc-kanMX6. The flexible Gs-linker was in- 
troduced in order to improve in vivo function of the epitope - 
tagged protein as suggested by Sabourin et al. (43). The con- 
struct was transformed using standard methods into strain 
FY4 (Y4015), a prototroph in the S288C background (44). 

The estradiol inducible MSN2 strains were constructed 
essentially as described in (45). First, MSN2-GFP was in- 
tegrated into the MSN2 locus of strain Y4098 (MATa 
msn2::URA3 leu2 AO ura3 AO) as described (35) to generate 
strain Y4105 MSW2 S288A ' S582A ' S620AS625AS633A ' S686A -GFP 
designated MSN2 6A , was constructed from the MSW2-GFP 
LEU2 plasmid B2802 using the QuikChange® II XL Site- 
directed mutagenesis kit (Agilent) and also integrated into 
strain Y4098 to generate strain Y4107. The Z 4 EV promoter 
was inserted upstream of MSN2-GFV or MSW2 6A -GFP lo- 
cus in these two strains by PCR amplification of the Z4EV 
promoter from plasmid pMNIO (45). The resulting strains, 
Y4108 and Y4109, were crossed with strain DBY12416 
(MATa leu2-F AC Ti-Z 4 E\ ybr034w::L£C/2) (45) and pro- 
totrophic segregants from each cross carrying the Z 4 EV 
transcription factor and the Z4EV-driven MSN2 locus were 
retained as Y4131 and Y4132. 



Growth protocols 

Glucose downshift. Yeast cells were grown at 30° C in Syn- 
thetic Complete (SC) + 2% glucose media and maintained 
in mid log phase for at least 24 h by dilution. When cul- 
tures reached a density of (3^1) x 10 6 cells/ml, cells were 
collected by filtration (Stericup-GP, 0.22-jjim filtering sys- 
tems, Millipore Corporation, Billerica, MA, USA), washed 
and resuspended in SC + 3% glycerol media and then sub- 
jected to sample preparation at the desired time points. 
Zero-minute samples were taken before this procedure. 

Estradiol induction. Cells were grown at 30° C in SD + 
2% glucose to a cell density of 4 x 10 6 cells/ml and estra- 
diol added to a final concentration of 10 |jlM. Samples (5 
ml) were harvested at specified times and ribonucleic acid 
(RNA) extracted as described below. 
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Table 1. 


Strains used in this study 




Strain 


Genotype 


Source 


Y2864 


MATol gall::HIS3 ade2-l canl-100 his3-ll,15 leu2-3,112 trpl-1 ura3-l 


(46) 


Y3513 


MATcl msn2 : : KanMX4 m.m4::KanMX4 gall::HIS3 ade2-l canl-100 


This study 




his3-ll,15 leu2-3,112 trpl-1 ura3-l (isogenic to Y2864) 




Y4105 


MATaleu2A0 ura3A0 MSN2-GFP 




Y4108 


MATaleu2A0 ura3A0 KanMX4-PzA E v-MSN2-GFP 


" 


Y4015 


MA 7a prototropli 


(44) 


Y4016 


AT^iraM5A'2-8xGly-13xMyc-KanMX4 (isogenic to Y4015) 


This study 


Y4131 


MA Taleu2A::P AC T] -Z 4 EV-NatMX KanMX4-Pz4 EV -MSW2-GFP 






ybr032w.:LEU2 




Y4132 


A#/<Ta/p«?A- , P^^T-<-7/FV-NatMX KanMX4-P7/cir-M < T/V7 6A -riFP 


„ 




ybr032w::L£CZ2 




Y4108 


MATa!eu2 A KanMX-Pz^v-MSWi-GFP 





Nucleosomal DNA sample preparation 

Mononucleosomal DNA from the wild-type parent strain 
and the msn2Amsn4A strain was isolated as previously 
described (3). Briefly, for each sample, a 650-ml culture 
was formaldehyde crosslinked, converted to spheroplasts, 
then resuspended in NP buffer, and micrococcal nuclease 
(Sigma-Aldrich, St. Louis, MO, USA) was added. DNA 
protein cross-links were reversed by incubation at 65°C for 
at least 4 h. DNA was then purified by PCR clean-up kit 
(Qiagen, Valencia, CA, USA), and the sample was analyzed 
by gel electrophoresis to ensure that the extent of diges- 
tion did not vary significantly from sample to sample. Data 
were obtained from single samples for each time point, with 
greater than 1.3 x 10 7 unique 2 x 100-bp reads per sample. 

mRNA sample preparation and transcript analysis 

Messenger RNA (mRNA) was isolated and hybridized to 
Agilent yeast microarrays as described in (46). Briefly, 5- 
ml cultures were collected on filters and snap frozen in liq- 
uid nitrogen. Total RNA was extracted using the Qiagen 
RNeasy Mini kit, including the additional DNase I diges- 
tion step. Chromosomal RNA (cRNA) for microarray hy- 
bridization was synthesized following the standard proto- 
col of the Agilent Low RNA Input Linear Amplification 
kit (Agilent Technologies). cRNA was extracted using the 
Qiagen RNeasy Mini kit and hybridized to Agilent Yeast 
Gene Expression Microarray (8 x 15K G4813A) slides and 
scanned at 5-|Jim resolution. Data were extracted using Ag- 
ilent Feature Extraction software version 9.5 with Linear 
Lowess dye normalization and no background subtraction 
and were submitted to the Princeton University Microarray 
database for storage and analysis. 

For estradiol induction experiments, time course fold 
change in transcript levels was fit to a Hill plot by optimiza- 
tion of n,fo, isTand F max for each gene for the equation f(t) — 
fo + V mAX -t n /(K 1 + t"). Delay times were determined by ex- 
trapolation of the derivative of this function aXj[t) — F max /2 
to the x-axis intercept. 

Chromatin preparation for chromatin immunoprecipitation 

Chromatin extract production was adapted from (47), with 
some modifications. Briefly, 45-ml yeast cultures prior to 
or post the glucose downshift procedure were crosslinked 
with formaldehyde (0.8% final concentration) for 10 min 



and quenched with glycine for 5 min. Cells were har- 
vested by centrifugation, 3000 X g, 4°C, 5 min, and 
washed with cold buffer (50-mM 4-(2-hydroxyethyl)-l- 
piperazineethanesulfonic acid (HEPES) pH 7.5, 140-mM 
NaCl), resuspended in 400-|xl cold ChIP lysis buffer (50- 
mM HEPES pH 7.5, 140-mM NaCl, 1-mM ethylenedi- 
aminetetraacetic acid (EDTA), 1% Triton, 0.1% sodium 
deoxycholate, 1-mM phenylmethylsulfonyl fluoride and a 
Roche complete protease inhibitor tablet) and snap frozen 
in liquid nitrogen. Samples were thawed in 37°C water bath, 
put on ice, and cold glass beads were added to 1 mm be- 
low meniscus. Cells were disrupted with a Fast Prep-24 (MP 
Biomedicals) bead beating system on setting 5.5 m/s 3 x 
40 s in a 4° C cold room. The resulting cell lysates were 
centrifuged at 20 800 x g, 4°C, 30 min. The supernatants 
were removed, and the pellets were resuspended in 100- 
|xl ChIP lysis buffer and placed in 120- (ulI Covaris tubes 
for sonication shearing. Chromatin was sheared to an av- 
erage fragment size of 350 bp using a Covaris E220 sys- 
tem. The sheared chromatin samples were transferred to an 
Eppendorf tube and sample volume adjusted to 200 |xl (by 
adding ChIP lysis buffer) and centrifuged at 10000 x g, 4°C, 
5 min. The pellets were the 'insoluble fraction' and the su- 
pernatants were transferred to a new Eppendorf tube and 
centrifuged again, 10 000 x g, 4°C, 15 min. The final super- 
natants were the chromatin extract used for ChIP. 



Chromatin immunoprecipitation 

For each ChIP, 2.5-jjlI anti-myc (Clontech, clone 9E10, 
cat#631206) or anti-Pol II C-terminal domain (Pol II 
8WG16 Monoclonal Antibody, Covance) antibody was 
added to 15-|xl resuspended protein G Dynabeads (Invit- 
rogen), coupled according to the Dynabeads manual and 
washed and resuspended in 233-|xl lysis buffer per sample. 
Sixty-seven microliter chromatin extract was incubated with 
the antibody-bound beads (total volume 300 |xl) with ro- 
tation for 4 h at room temperature (RT). The beads were 
then collected with the magnet and washed (resuspended 
and nutated 4 min, RT) with 0.5-ml sodium dodecyl sul- 
phate (SDS) buffer (50-mM HEPES pH 7.5, 140-mM NaCl, 
1-mM EDTA and 0.025% SDS). Beads were subsequently 
washed with 0.5-ml high salt buffer (50-mM HEPES pH 
7.5, 500-mM NaCl, 1-mM EDTA), followed by 0.5-ml tris- 
lithium (TL) buffer (20-mM Tris-Cl pH 7.5, 140-mM NaCl, 
250-mM LiCl, 1-mM EDTA), followed by two washes in 
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0.5-ml tris-EDTA (TE) buffer (20-mM Tris-Cl pH 7.5, 0.1- 
mM EDTA). Washed beads were resuspended for elution 
in 72-jxl TE + 1% SDS buffer (20-mM Tris-Cl pH 7.5, 0.1- 
mM EDTA, 1% SDS), vortexed and heated in a 65°C water 
bath, 2 min. The beads were vortexed well again and super- 
natants were taken from the beads. Twenty-five microliter 
was used for western blots and 45 |jl1 was taken to reverse 
crosslink at 65°C, 14 h. 



Reverse crosslinking and purification of DNA 

Input DNA (2-(xl chromatin extract (input DNA) and 118- 
ixl TE + 1% SDS buffer) and ChIP DNA (45-julI ChIP eluate 
+ 75-|xl TE + 1% SDS buffer) were incubated at 65°C, 14 h 
for reverse crosslinking. Reverse crosslinked samples were 
purified on Qiagen PCR purification columns, eluted in 2 
x 35-|xl Qiagen Elution buffer and kept frozen until library 
construction. 



Library construction 

ChlP-DNA was amplified using the LM-PCR method de- 
scribed in Agilent Yeast ChlP-on-chip analysis protocol 
version 9.2, May 2007 and subjected to the Illumina TruSeq 
paired-end sequencing protocol. 



Sequence analysis 

Paired-end sequences were mapped to the Saccharomyces 
cerevisiae reference genome s288c Saccharomyces Genome 
Database (SGD) version r64-l-l .20110203, using Bowtie 
for Illumina (version 1.1.2) with seed length 22 and a maxi- 
mum permitted total of quality scores of 70 at mismatched 
read positions, also allowing a maximum of two mismatches 
in the seed. Twenty samples each for the Msn2 ChIP at 0 and 
20 min were combined to yield 2 x 10 5 total reads for each 
time point. Alignments that mapped to more than one posi- 
tion on the reference genome were randomly distributed be- 
tween the reportable alignments. To eliminate PCR amplifi- 
cation artifacts, precise duplicates of paired-end ChIP read 
alignments mapping to a genomic position were excluded 
from analysis. Resulting sequence positions were then sub- 
jected to further analysis in MATLAB. Occupancy at each 
base pair position across the genome for both nucleosomes 
and ChIP profiles was determined by summing the total 
number of unique sequence reads at that position and then 
normalizing the summed values such that the average occu- 
pancy per bp for each experiment equals 1 over each chro- 
mosome. Peaks of Msn2 binding were identified either as 
those with a maximal peak intensity 6-fold above the av- 
erage binding over the chromosome in which it is located 
or as those with z-score greater than 5 for the integrated 
area of binding in the 250-bp region around a binding maxi- 
mum. Visualizations were performed using MATLAB stan- 
dard bioinformatics methods. The positions of STRE ele- 
ments were obtained from SGD (http://www.yeastgenome. 
org/cgi-bin/PATMATCH/nph-patmatch). Functional anal- 
ysis of groups of genes was performed using the Gene On- 
tology Term Finder from SGD. 



Accession numbers and deposition of microarray data 

Read data for the ChlP-Seq and MNase-Seq experiments 
are publically available at NCBI SRA with the acces- 
sion number SRP033438. Microarray data are publicly 
available at http://puma.princeton.edu/cgi-bin/publication/ 
viewPublication.pl?pub_no=559 and as a processed spread- 
sheet in Supplementary Table S2. 

RESULTS 

Msn2 binds to a limited number of sites in vivo 

To explore the relation between transcription factor bind- 
ing, transcriptional changes and nucleosome repositioning, 
we determined the global binding pattern of Msn2 by chro- 
matin immunoprecipitation and DNA sequencing of the 
precipitated fragments (ChlP-Seq) prior to and 20 min after 
transition of cells from growth on glucose to growth on glyc- 
erol, a condition that induces the ESR. We performed ChlP- 
Seq using anti-Myc antibodies on a strain in which MSN2 
was replaced with MSN2 tagged with 13 copies of the Myc 
epitope attached to the carboxy terminus of the protein and 
expressed under its own promoter. The Myc-tagged version 
of the protein showed normal nuclear localization and tran- 
scriptional activation in response to both hydrogen peroxide 
and glucose downshift conditions (Elfving et al. 2014, sub- 
mitted). We obtained 3-4-fold average sequence coverage 
over the entire genome for both time points and 180 reads 
over the most abundant unique binding site at the 20 min 
time point. To assess the interplay of nucleosome remodel- 
ing and Msn2 binding, we concurrently mapped genome- 
wide nucleosome positions prior to and 20 min after the 
glucose-to-glycerol switch in an MSN2 MSN4 strain and 
in an isogenic msn2 msn4 strain by sequencing size-selected 
DNA fragments following micrococcal nuclease treatment 
of cross-linked chromatin. We obtained > 100-fold sequence 
coverage of the entire genome for both strains at each time 
point. 

ChlP-Seq identified few Msn2 binding sites prior to the 
carbon source downshift and a large number after the 
downshift. We computationally identified sites of Msn2 
binding as described in the Materials and Methods section. 
The positions of the major Msn2 binding sites are shown in 
Figure 1. We hand annotated each of the peaks to identify 
the genomic features associated with each site. This process 
yielded 273 distinct and robust peaks of bound Msn2, dis- 
tributed over 269 genes, 20 min after the glucose downshift. 
The positions of these sites, the associated gene or genomic 
feature and the relative abundance of Msn2 at these sites 
prior to and after the glucose downshift are listed in Sup- 
plementary Table SI. The majority (192) of those sites cor- 
responded to promoter regions with the remaining peaks 
mapping solely to gene coding sequences (36 peaks) or to 
Ty elements (40 peaks), with the latter displaying a very 
characteristic pattern of Msn2 distribution based on the 
manner in which sequence reads were apportioned to re- 
peat sequences. Because Ty elements are not readily distin- 
guished by sequence, we could not determine whether all Ty 
elements bind Msn2 at equal levels or whether some have 
greater affinity than others. Most of the coding sequences 
registering significant Msn2 binding were expressed at high 
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Figure 1. Msn2 binding sites. The relative positions of the 268 Msn2 binding sites determined in this study are indicated by short vertical lines above 
each chromosome (horizontal black lines). The vertical lines below the chromosomes denote the 212 most robust binding sites identified in (41) following 
exposure of cells to hydrogen peroxide. The binding sites identified in the current study are denoted by color specified in the legend as residing solely in the 
coding region of a gene (coding region), over a transposable element (Ty) or in the promoter of a gene induced, repressed or unaffected (neutral) by Msn2. 
See Supplementary Table SI for a detailed description of each site. 



levels, as measured by PolII occupancy (Supplementary Ta- 
ble SI ), consistent with the growing appreciation that highly 
expressed genes are retrieved inadvertently as artifacts of 
the ChIP protocol (48). In fact, more than half of the 50 
most highly expressed genes were recovered in the Msn2 
ChIP experiment (P < 10~ 100 ). Moreover, most coding re- 
gions to which Msn2 bound lack an STRE, the binding mo- 
tif for Msn2. In sum, our experiments identified a number 
of sites for Msn2 binding following glucose downshift, with 
a majority of the robust sites mapping to promoter regions 
and a smaller number mapping to Ty elements and to cod- 
ing domains of highly expressed genes. 

Previous in vivo and in vitro studies have identified a 
canonical binding site for Msn2, the STRE with a sequence 
RGGGG (26,37,49-51). As shown in Figure 2A, Msn2 
binding is significantly enriched around STREs. Moreover, 
the probability of an STRE residing close to a peak of Msn2 
binding is significantly greater than that expected for a ran- 
dom distribution of STREs relative to binding sites (Fig- 
ure 2B). However, while approximately 11,500 STRE se- 
quences are present in yeast genome, with 3160 lying in 
the promoters of almost 2000 genes, Msn2 fails to bind 
most of these sites. Reasons for this selective binding to 



only a subset of STREs are addressed below. The proba- 
bility of Msn2 binding to a promoter was correlated with 
the number of STREs within that promoter (Supplemen- 
tary Figure SI). This may result from cooperative binding 
of Msn2 to adjacent STREs, particularly since the fraction 
of multiple STRE-containing promoters that are bound by 
Msn2 is higher than that predicted assuming independent 
interaction of Msn2 with each STRE within a promoter 
(Supplementary Figure SI). We also observed low but de- 
tectable levels of Msn2 binding over STREs prior to the 
nutrient downshift. This is consistent with the observation 
that, while Msn2 resides predominantly in the cytoplasm in 
the absence of stress, Msn2 is not completely excluded from 
the nucleus under those conditions (35). Moreover, even un- 
der robust growth conditions, Msn2 exhibits random nu- 
clear bursting, such that Msn2 congregates in the nucleus 
for a short period of time in a small number of individual 
cells. Thus, the low level of Msn2 binding we observe on av- 
erage in unstressed cells may represent robust Msn2 binding 
in a very small number of cells in the population. 

We compared the binding sites identified in our experi- 
ment with those observed by Huebert et al. following treat- 
ment of cells with hydrogen peroxide (41). This study re- 
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Figure 2. Msn2 binds near STREs. (A) The average Msn2 occupancy at a particular site is plotted as a function of distance of that site from each of the 
3150 promoter STREs (RGGGG motifs), before and 20 min after a glucose-to-glycerol downshift. (B) The likelihood of an STRE residing at a particular 
site is plotted as a function of distance of that site from each of the peaks of Msn2 binding identified in this study and listed in Supplementary Table SI. The 
values are calculated as the fraction of Msn2 peaks that contain an STRE at the indicated position, averaged over a 20 bp moving window. By comparison, 
the probability of finding STRE at a random site is less than 9 x 10 -4 . 



ported 1289 sites that were statistically enriched in the 
treated samples. Plotting the 212 sites with the highest re- 
ported peak values from the Huebert et al. experiment rel- 
ative to the positions of the most robust sites from our ex- 
periment indicates a significant overlap in the two sets of 
data (Figure 1). In fact, the peaks of 113 sites (42%) iden- 
tified in our experiment lie within 500 bp of a peak among 
the 212 most abundant sites identified by Huebert et al. (P 
< 10~ 250 ). Moreover, 76% of the peaks identified in our ex- 
periment coincide with the location of one of the significant 
sites identified by Huebert et al. Thus, the vast majority of 
de novo Msn2 binding sites are recapitulated under different 
stress conditions, consistent with a relative uniform tran- 
scriptional response of cells to stress, independent of the 
nature of the initiating stress. 

Stress-specific binding of Msn2 to promoters reflected 
both common and condition-specific responses to stress. 
While most of the Msn2-binding sites identified upon nu- 
trient downshift overlap those identified following oxida- 
tive stress, we noted that a number of stress-induced Msn2 
binding sites identified after oxidative stress were absent 
upon nutrient downshift, and vice versa. Those genes in 
whose promoters Msn2 bound following oxidative stress 
but not after glucose downshift were enriched for telom- 
eric sequences and for those involved specifically in response 
to treatment with oxidizing agents (13/69; P — 10~ 4 ). Sim- 
ilarly, those genes in whose promoters Msn2 bound fol- 
lowing glucose downshift but not in response to oxidative 
agents were highly enriched in those involved in carbohy- 
drate metabolism (P=4 x 10" 5 ). Thus, in addition to bind- 
ing to promoters of genes involved in a common stress re- 
sponse, Msn2 binds specifically to a subset of genes asso- 
ciated with response to the particular initiating stress. We 
discuss below possible mechanistic bases for this discrimi- 
nation. 

Msn2 mediates both gene activation and gene repression 

We assessed the effect of Msn2 on transcriptional repro- 
gramming under nutrient downshift conditions in several 
ways. First, we determined the level of Pol II associated 



with all coding regions prior to and 20 min after glucose 
to glycerol transition by performing ChlP-Seq using anti- 
bodies targeting the C-terminal repeat domain of the Rpbl 
subunit of Pol II. In addition, we examined global transcript 
levels by microarray analysis of both an MSN2 MSN4 and 
an msn2 msn4 strain pre- and post-transition. These data 
provide information on the extent to which transcriptional 
changes that occur following the carbon source downshift 
are dependent on the Msn transcription factors. Finally, to 
identify those transcriptional changes that occur specifically 
in response to activation of Msn2, we measured global tran- 
scriptional changes following ectopic induction of Msn2. 
This was accomplished using a hybrid zinc-finger transcrip- 
tion factor recently described, termed Z4EV (45). Specifi- 
cally, we used a strain in which MSN2 expression was driven 
by a modified GAL1 promoter in which four repeats of 
a zinc-finger DNA binding protein recognition sequence 
replaced the Gal4 binding sites. The strain expressed the 
Z 4 EV fusion protein, consisting of the zinc -finger binding 
protein targeting the modified GAL1 promoter as well as the 
VP 16 activation domain and an estrogen receptor fragment 
encompassing the Hsp90 and estrogen binding domains. 
Z 4 EV resides in the cytoplasm in the absence of estradiol, 
due to binding Hsp90, and rapidly dissociates from Hsp90 
following estradiol addition, resulting in relocation to the 
nucleus and induction of transcription of genes with a up- 
stream activator sequence (UAS) containing the specific 
zinc-finger binding motifs. Since this sequence does not nor- 
mally exist in the yeast genome, the only gene induced by 
estradiol treatment is the one engineered to be linked to that 
sequence. Accordingly, estradiol induces production only of 
Msn2 in this strain. While the induced Msn2 resides pre- 
dominantly in the cytoplasm upon estradiol treatment, suf- 
ficient nuclear localization occurs, through partitioning in 
all cells and bursting in some cells, to allow Msn2 to ex- 
ert its transcriptional effects (see below). We conclude that 
any gene rapidly induced or repressed by estradiol treatment 
must be under the direct regulation of Msn2. 

Using the transcript data from the experiments described 
above, we could identify those genes whose expression 
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was directly affected by Msn2 in response to the glucose 
downshift. The strong correlation between transcript level 
changes and the changes in Pol II occupancy over the corre- 
sponding coding regions following nutrient downshift con- 
firmed that transcript level changes were a consequence 
of changes in transcriptional activation rather than post- 
transcriptional processes (Supplementary Figure S2). We 
identified those genes activated by Msn2 as those that 
showed increased transcript levels upon estradiol treatment 
of the Z4EV strain described above as well as diminished in- 
duction, or more substantial repression, of transcript levels 
in the msn2 msn4 strain versus the MSN2 MSN4 strain af- 
ter the glucose downshift. Similarly, we identified genes re- 
pressed by Msn2 as those whose transcript levels fell upon 
estradiol treatment of the Z4EV strain and exhibited higher 
transcript levels in the msn2 msn4 strain versus the MSN2 
MSN4 strain after glucose downshift. These independent 
measures of sufficiency and necessity of Msn2 activity on 
gene expression were reasonably consistent (Supplementary 
Table SI). Moreover, approximately two-thirds of genes 
exhibiting Msn2-dependent regulation by the above crite- 
ria showed Msn2 promoter binding in response to a glu- 
cose downshift, consistent with the hypothesis that Msn2 
binding directly affected expression of the corresponding 
gene and that most genes whose transcription modula- 
tion are Msn2-dependent are directly regulated by Msn2 
binding. Finally, a significant number of genes (86/192) to 
whose promoter Msn2 bound after the glucose downshift 
showed Msn2-dependent transcriptional activation, consis- 
tent with Msn2's reported role as a transcriptional activa- 
tor. These were enriched in genes involved in energy reserve 
metabolism (P — 3 x 10 -11 ), oxidation-reduction processes 
(P = 5 x 10" 10 ) and glycogen (P = 7 x 1(T 10 ) and trehalose 
(P — 4 x 10 -7 ) metabolism. However, a significant num- 
ber of genes (44/192) to whose promoter Msn2 bound ex- 
hibited Msn2-dependent transcriptional repression follow- 
ing glucose downshift or during induction of Msn2. These 
were enriched in genes involved in glucose catabolism (P — 
10~ 4 ). The remaining genes to which Msn2 bound were ei- 
ther Ty elements or coding regions noted above or showed 
no Msn2-dependent change in expression. These results in- 
dicate that Msn2 functions both as a transcriptional acti- 
vator and a transcriptional repressor. The basis of this dual 
activity is discussed below. 

Msn2 elicits different patterns of gene regulation kinetics 

Evaluation of the transcriptional consequences of activat- 
ing Msn2 using the Z4EV system revealed several unex- 
pected aspects of Msn2 regulation. First, a number of genes 
significantly changed expression following Z4EV induction 
of Msn2 but did not exhibit significant Msn2 binding in 
the ChlP-Seq analysis or display Msn2-dependent tran- 
scriptional changes following the glucose downshift. The 
induced genes in this set generally contained one or more 
STREs in their upstream intergenic regions. For instance, 
85 of the 100 most induced genes following estradiol treat- 
ment of the Z4EV strain contained one or more upstream 
STREs, although only 40 of these showed significant bind- 
ing of Msn2 by ChlP-Seq in the glucose downshift experi- 
ment. This suggests a hierarchy of STRE binding affinities 



such that lower affinity sites are bound only when Msn2 is 
expressed at higher levels but that such sites can mediate ac- 
tivation under that condition. A number of these promoters 
exhibited Msn2-dependent nucleosome remodeling (see be- 
low) during the glucose downshift, indicating that even low 
level or transient Msn2 binding can affect local chromatin 
structure. Those genes repressed following Msn2 induction 
by Z 4 EV are highly enriched for ribosome biogenesis genes 
(P = 10~ 24 ) but only -10% show Msn2 binding by ChlP- 
Seq and most (60%) do not contain STREs in their promot- 
ers. This suggests that much of the repression is an indirect 
effect of Msn2 induction. 

The Z4EV induction data also revealed unexpected diver- 
sity in the kinetics of activation of Msn2 responsive genes. 
Genes induced upon activation of Msn2 in the Z4EV strain 
begin to accumulate transcripts following estradiol addition 
but only after a delay. As illustration, the induction kinetics 
of several genes are shown in Figure 3A-D. The duration of 
this delay showed a broad distribution among the induced 
genes, with the majority of genes initiating transcript accu- 
mulation 35-45 min after estradiol addition while a signif- 
icant minority of the genes initiated accumulation greater 
than 45 min following addition (Figure 31 and J). Since the 
Z4EV cells in which Msn2 was induced by estradiol were un- 
stressed, Msn2 remained predominantly in the cytoplasm, 
with only occasional random bursts of nuclear occupancy in 
individual cells (Supplementary Movie SI and Supplemen- 
tary Figure S3A). To test whether this limited nuclear occu- 
pancy contributed to the variable lag in transcript accumu- 
lation, we performed the transcription study using the Z 4 EV 
strain to drive expression of a mutant Msn2, designated 
Msn2 6A , in which all six PKA phosphorylation sites in the 
protein were converted to alanines. Since PKA-mediated 
phosphorylation is responsible for retaining Msn2 in the 
cytoplasm and Msn2 6A could not be phosphorylated by 
PKA, Msn2 6A induced by estradiol treatment immediately 
entered the nucleus in all cells (Supplementary Movie S2 
and Supplementary Figure S3B). Notably, the duration of 
the delay in transcript accumulation was reduced for most 
genes to 15-25 min (Figure 3E-J), without a significant 
change in the rate of transcript accumulation once induc- 
tion began (Figure 3K). Repressed genes exhibited a similar 
pattern of kinetics (data not shown). Thus, Msn2 responsive 
genes fall into roughly two categories based on their activa- 
tion kinetics in conditions of limited Msn2 nuclear occu- 
pancy: rapid responders and slow responders. Notably, this 
dichotomy is lost under conditions of high Msn2 nuclear 
occupancy, ruling out the possibility that the two classes 
simply reflect direct versus indirect targets of Msn2 regu- 
lation. Rather, Msn2-regulated genes exhibit distinctly dif- 
ferent patterns of response to wild-type Msn2 activation. 

Msn2 can promote repositioning of nucleosomes following nu- 
tritional stress 

Consistent with our previous observations (3), we find that 
nutrient downshift results in relatively limited nucleosome 
repositioning, in spite of the significant transcriptional re- 
programming (Supplementary Figure S4). However, by de- 
termining global nucleosome positions before and after glu- 
cose downshift in both an MSN2 MSN4 and an msn2 msn4 
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Figure 3. Different genes exhibit different induction kinetics in response to Msn2. (A-D) Induction kinetics for four Msn2-inducible genes following 
estradiol addition to strains containing wild-type MSN2 under control of the hybrid Z4EV transcription factor. Blue circles: fold increase in gene transcript 
levels relative to that immediately prior to estradiol addition; solid red line: best fit of the data for each gene to the function f[t) =fo+ K max ■ t"/(K" + 
/"); dashed line: tangent line to the curve atf(t) = K max /2, whose extrapolation to the x-axis provides the measure of the time delay in response to Msn2 
induction. (E-H) Induction kinetics for the genes shown in (A-D) following estradiol addition to strains containing MSN2 6A under control of the hybrid 
Z4EV transcription factor. (I) Histogram of time delay values for the 96 genes induced more than 2-fold in both the MSN2 wild-type and MSN2 6A strains 
and whose induction values are reasonably fit by the Hill curve. Blue line: delay values in the MSN2 wild-type strain; red line: delay values in the MSN2 6A 
mutant strain. (J) Scatter plot of the delay time for each gene in I in the MSN2 6A strain relative to that in the MSN2 wild-type strain. (K) Histogram of 
rates of induction, i.e. the slope of the tangent line to the fitted curve at f(t) = F max /2, for the 96 genes in (I). Blue line: delay values in the MSN2 wild-type 
strain; red line: delay values in the MSN2 6A mutant strain. 
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Figure 4. Msn2 promotes nucleosome repositioning over gene promoters. Shown are binding profiles of Msn2 (row 1) and PolII (row 2) as well as the 
nucleosome profiles (rows 3 and 4) over four different Msn2-regulated genes and the surrounding regions along the genome. Each panel shows a binding 
profile before the glucose-to-glycerol switch (blue lines) and 20 min after the transition (red lines). The nucleosome profiles were obtained for both an 
MSN2 MSN4 (row 3) and an msnl msn4 strain (row 4). Two of the genes, PGM2 and HSP30, are induced by Msn2 following the nutrient downshift, while 
the other two, ERG11 and ZRT1, are repressed. In order to be able to compare the occupancy profiles obtained from different experiments with different 
sequencing coverage, we normalized all profiles such that the average occupancy across each chromosome in each experiment equals 1. 
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Figure 5. Msn2 reduces the average nucleosome occupancy near its binding sites following nutrient downshift. The distribution of nucleosome occupancy 
around STREs (A) and Msn2 binding sites (B), as defined in the legend to Figure 2, are shown before (blue line) and after (green line) the glucose-to-glycerol 
downshift in wild-type cells and before (red line) and after (cyan line) the glucose-to-glycerol downshift in msn2 msn4 mutant cells. 



strain, we were able to identify promoters in which expres- 
sion change was associated with nucleosome repositioning 
and determine whether that repositioning was dependent 
on Msn2 and/or Msn4. In this manner, we identified a num- 
ber of genes in which transcriptional activation, Msn2 bind- 
ing and nucleosome depletion from the promoter were co- 
incident and in which nucleosome depletion was dependent 
on Msn2 and/or Msn4. Two examples of genes, PGM2 and 
HSP30, that show this pattern are profiled in Figure 4A 
and C. Similarly, we identified several genes at which tran- 



scriptional repression, Msn2 binding and nucleosome ac- 
quisition are coincident and in which nucleosome acquisi- 
tion is dependent on Msn2 (Figure 4B and D). In sum, 72% 
(62/86) of induced genes regulated by Msn2 show nucleo- 
some depletion and 60% of those show complete or partial 
dependence of the nucleosome depletion on Msn2 and/or 
Msn4. Similarly, 50% (22/44) of repressed genes regulated 
by Msn2 show nucleosome remodeling, primarily nucleo- 
some acquisition, and 45% of those show dependence of 
nucleosome remodeling on Msn2 and/or Msn4. From these 
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Figure 6. Msn2 binds to STREs in nucleosome free promoters. All STREs residing upstream of gene transcription start sites were organized by ^r-means 
clustering (k = 10) on the basis of the nucleosome occupancy profiles over the 1 kb regions centered on STRE at 20 min following the glucose-to-glycerol 
nutrient downshift. The clustered nucleosome occupancy profiles are presented as a heat map (A). A heat map of the binding profiles of Msn2 at 20 min 
following the nutrient downshift is shown in (B), with the same gene order as in (A). The average nucleosome and Msn2 occupancy in the lkb regions 
centered on STRE elements for each cluster, before and 20 min after the switch, are shown in (C). 
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Figure 7. Msn2 and STRE distribution relative to the well-positioned nu- 
cleosomes. We selected well-positioned nucleosome, with dyad positions 
characterized by at least five sequencing reads. The blue line shows the av- 
erage Msn2 occupancy relative to the dyad locations of the well-positioned 
nucleosomes. The green line shows the probability of finding a nearby 
STRE element. 



by other agents but plays an active role in restructuring the 
NDR during transcriptional reprogramming. 

The Msn2/Msn4-dependent reorganization of nucleo- 
somes noted above could be observed on a global scale. We 
calculated the nucleosome occupancy around all promoter 
STRE sites before and after the nutrient downshift in both 
the wild-type and msn2 msn4 strains. As evident from the 
plot of the average nucleosome occupancy as a function of 
distance from every STRE under these four conditions, nu- 
trient downshift results in reduction of nucleosome density 
over these sites (Figure 5A). This is consistent with an over- 
all increase in expression of stress-responsive genes follow- 
ing a downshift, and observed correlation between reduced 
occupancy of nucleosomes in promoters and increased tran- 
scription. However, the reduction in nucleosome density of 
STREs was approximately twice as large in the wild-type 
strain compared to the msn2 msn4 strain. We observed sim- 
ilar results by examining nucleosome density around sites 
of Msn2 binding (Figure 5B). These observations demon- 
strate that Msn proteins play an active role in reducing nu- 
cleosome occupancy during transcriptional activation. 



studies, we conclude that a significant function of Msn2 is to 
expose promoter regions during gene activation and to oc- 
clude promoter regions during gene repression. Thus, Msn2 
is not simply a passive respondent to remodeling promoted 



Nucleosomes restrict access of Msn2 to STREs 

While our data described above demonstrate that Msn2 
binding can alter adjacent nucleosome occupancy, we find 
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Figure 8. Possible mechanisms for stress-specific binding of Msn2 to different sets of genes. (A) Different stresses could result in distinct modifications (red 
versus green dots) of Msn2 (blue oval), which could alter the binding specificity or the nuclear occupancy dynamics and which would lead to interaction 
with different but overlapping sets of stress-responsive genes. (B) Different stresses could activate Msn2 (blue oval) as well as a stress-specific transcription 
factor (red oval for nutrient stress and green oval for oxidative stress). Those genes with STREs lying in nucleosome-free domains would bind Msn2 under 
either condition. However, binding of the stress-specific transcription factor could partially unwrap adjacent nucleosomes (beige ovals) to reveal additional, 
previously inaccessible STREs to which Msn2 could bind. 



that nucleosome position significantly restricts accessibil- 
ity of Msn2 to its canonical binding sites. First, while most 
STREs in the genome reside in coding regions, Msn2 did 
not bind to any of these sites, likely because they are oc- 
cluded by well-positioned nucleosomes. Second, we used k- 
means clustering to categorize patterns of nucleosome oc- 
cupancy around the 3150 transcription start site proximal 
STREs prior to the nutrient downshift (Figure 6A and C). 
These patterns range from cases in which nucleosomes sit 
directly over the STRE to cases in which the STRE resides 
in an NDR centered on the STRE to cases in which the 
STRE resides in an extended NDR. We then determined 
to which of these STREs Msn2 bound following the nu- 
trient downshift. As evident in Figure 6B and C, Msn2 
bound almost exclusively only to those STREs residing in 
extended NDRs. Moreover, the extent of Msn2 binding was 
essentially inversely proportional to the nucleosome den- 
sity around the STRE. Surprisingly, even those STREs lying 
in open chromatin domains of limited dimension were not 
substrates for Msn2 binding (cluster 2, for example). Thus, 
at the global scale, nucleosome occupancy restricts Msn2 
binding. 

We also find a local effect of nucleosome occupancy on 
Msn2 binding. In Figure 7 we plot the density of STRE ele- 
ments within the footprint of well-positioned nucleosomes 
across the genome. Superimposed on that plot is the aver- 
age binding of Msn2 to STREs at the indicated position 
following nutrient downshift. Quite evident from this dia- 
gram, STREs are enriched under the nucleosome umbrella, 
perhaps reflecting the guanosine/cytosine bias in nucleo- 
somal positioning sequences. In contrast, Msn2 binding is 
uniformly low within the 90-bp inner core of the positioned 
nucleosomes, significantly higher in the nucleosome adja- 
cent region and proportionately increasing with increasing 
distance from the inner core of the nucleosome footprint. 
These results suggest that well-positioned nucleosomes pre- 
vent access of Msn2 to its cognate binding sites when those 
sites reside under the core of the nucleosome. However, 
Msn2 can gain access to its cognate sites that lie under the 
edges of the positioned nucleosomes. As discussed below, 
this observation is consistent with dynamic partial unwrap- 



ping of DNA bound to nucleosomes and competition be- 
tween nucleosome binding and transcription factor binding 
to specific target sequences. 

DISCUSSION 

Stress-specific and stress-non-specific Msn2 binding 

We have mapped Msn2 binding sites across the genome fol- 
lowing a nutrient downshift, which elicits the ESR. Some 
of the binding sites correspond to structural features such 
as transposable elements or transcribed regions of highly 
expressed genes. The former are repeated sequences, which 
prevents assignment of binding to specific elements, so we 
cannot determine from our data whether all Ty elements 
or only a subset binds Msn2. However, binding to Ty ele- 
ments is consistent with stress activation of Ty transposi- 
tion (52,53) and suggests that Msn2 may play a role in stim- 
ulating such transposition. The latter structural feature- 
-transcribed coding regions — may be an artifact of the 
ChlP-Seq protocol, since antibodies to many unrelated pro- 
teins have been reported to retrieve these same sequences 
(48,54,55). We also find that anti-Myc antibodies to mul- 
tiple different tagged proteins in this experiment retrieved 
these sequences (data not shown). The majority of the bind- 
ing sites correspond to promoters of transcribed genes. 
Moreover, most of these binding sites lie within several hun- 
dred base pairs of one or more STREs, previously defined 
as Msn2 binding sites. We conclude that most Msn2 bind- 
ing is targeted to specified Msn2 binding motifs in promoter 
regions of genes. 

Comparing our data with those previously obtained for 
Msn2 binding following peroxide treatment revealed a com- 
mon core of Msn2-bound genes as well as condition-specific 
binding. Since we examined binding at only one time point 
following nutrient shift, we may not have captured all the 
nutrient-dependent binding sites. Nonetheless, the bind- 
ing pattern mirrors the transcriptional changes associated 
with different stresses, in which a core stress response is 
augmented by activation/repression of genes targeted to 
specific stresses (20,21). Consistent with this interpreta- 
tion, those genes bound by Msn2 uniquely following per- 
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oxide treatment are enriched for those associated with re- 
mediation of oxidative stress while those uniquely bound 
following glucose downshift are involved in carbohydrate 
metabolism. These results suggest that the stress-specific 
transcriptional response is mediated at least in part by 
Msn2. 

How might Msn2 binding be responsive to specific stim- 
uli? One possibility is that post-translational modification 
of Msn2 dictated by different stress signaling pathways 
modifies the binding specificity of the protein (Figure 8A). 
In this context, different stresses yield different patterns of 
Msn2 nuclear entry and exit, patterns that may well be per- 
ceived by different promoters in different ways (34-36). We 
found in this study that Msn2 responsive genes had quite 
divergent delay times in response to Msn2 activation when 
Msn2 was largely cytoplasmic and exhibited only random 
bursts of nuclear occupancy. This difference was essentially 
eliminated when Msn2 resided predominately in the nu- 
cleus. This divergence in response times may reflect the dis- 
tinction recently described by Hansen and O'Shea (36), who 
noted that some genes respond to both sustained and pul- 
satile Msn2 nuclear localization (fast promoters) while oth- 
ers responded only to sustained Msn2 localization (slow 
promoters). In fact, two of the three slow-responding genes 
defined by Hansen and O'Shea, ALD3 and TKL2, were 
among the slowest responding genes to wild-type Msn2 in- 
duction in our study, while all four of the fast-responding 
genes from Hansen and O'Shea, DDR2, DCS2, RTN2 and 
HXK1, initiated mRNA accumulation significantly earlier, 
along with the majority of induced genes (Figure 3A-D). 
The third slow-responding gene characterized by Hansen 
and O'Shea was not captured adequately for analysis in our 
microarray study. We could not identify a correlation be- 
tween any aspect of nucleosome remodeling with response 
delay times: both fast- and slow-responding genes exhibited 
Msn2/4-dependent clearance of nucleosomes from their 
NDRs, for instance. However, further analysis of Msn2 
binding, response kinetics and nucleosome remodeling may 
prove informative. 

Another possible scenario to explain stress-specific Msn2 
promoter binding is that transcription factors responsive 
to specific stress signals, e.g. Yapl for oxidative stress or 
the Hapl-4 complex for glucose downshift (56-59), bind to 
promoters in a stress-specific manner and stimulate chro- 
matin clearance at those promoters, allowing access of 
Msn2 to STRE sites that would otherwise be occluded (Fig- 
ure 8B). This may represent an example of cooperative bind- 
ing of transcription factors by sequential unwrapping of 
DNA from nucleosomes, such that binding of a transcrip- 
tion factor to its cognate site near the periphery of a po- 
sitioned nucleosome provides access to a binding site for a 
second transcription factor that would be otherwise buried 
under the interior of the nucleosome (60,61). Consistent 
with this model, we note a significant overlap (P < 10~ 24 ) 
between Msn2 binding and promoters at which Floer et al. 
mapped RSC-associated and partially unwrapped nucleo- 
somes (62). Moreover, we find that more than 70% of the 
promoters to which Msn2 binds and activates transcription 
undergo nucleosome remodeling and for 40% of those the 
remodeling is independent of Msn2. Thus, other transcrip- 



tion factors may well clear the space to allow Msn2 binding 
and that clearance may well be stress specific. 

Msn2 promotes both transcriptional activation and transcrip- 
tional repression 

We find that Msn2 binding stimulates both transcriptional 
activation and transcriptional repression. The capacity of 
Msn2 to promote transcriptional activation is well docu- 
mented and consistent with the structural features of the 
protein (26,63). The activity as a repressor is less well docu- 
mented. Our data demonstrate that repression is not an in- 
direct effect, as might result from transcriptional activation 
of a repressor protein or inhibition of growth. Rather, Msn2 
binds to promoters of repressed genes and in some cases is 
responsible for recruitment of nucleosomes into the NDR. 
How Msn2 binding results in activation in some cases and 
repression in others is certainly not clear but may involve 
the type of combinatorial interaction with transcriptional 
modulators as mentioned above. Moreover, in a companion 
paper (Elfving et al. 2014, submitted), we show that Msn2 
recruits mediator complex, most often to promote recruit- 
ment and activation of Pol II but occasionally to reposition 
Mediator to a non-productive position in the promoter . 
Thus, the same basic activity can function both in activa- 
tion and repression. Finally, since stress is associated with 
at least a transient cessation of growth (58,64,65), we were 
interested in understanding whether Msn2 might repress 
genes whose expression is necessary for growth. In fact, a 
significant component of the ESR consists of repression of 
genes that promote growth, such as ribosomal protein and 
ribosome biogenesis genes (21). We do find that many of 
the genes repressed upon activation of Msn2 are highly en- 
riched for those involved in ribosome biogenesis. However, 
few of these genes are bound by Msn2, at least under condi- 
tions of nutrient downshift. Rather, we observed that Msn2 
activates transcription of DOT1, which encodes a repressor 
of ribosome biogenesis genes (see Supplementary Tables SI 
and S2). Moreover, we find that Msn2 binds to and acti- 
vates transcription of XBP1, which encodes a repressor of 
a number of genes required for cell cycle progression (66). 
Accordingly, Msn2, while a primary purveyor of the ESR, 
may indirectly repress the growth-associated genes encom- 
passed in the ESR. 

A complex interplay between Msn2 binding and nucleosome 
occupancy 

More than 10,000 canonical Msn2 recognition sites reside in 
the yeast genome and yet only a small fraction of these serve 
as binding sites for Msn2 in vivo. Comparing nucleosome 
occupancy to subsequent Msn2 binding, we see that those 
STREs that fail to serve as binding sites generally lie in re- 
gions of well-ordered nucleosomes. Moreover, those STREs 
lying under the core of a well-positioned nucleosome show 
diminished binding of Msn2, relative to sites outside nucle- 
osomes or at the edges of positioned nucleosomes. Thus, po- 
sitioned nucleosomes serve to restrict Msn2 binding. More- 
over, the gradient of Msn2 binding as a function of the dis- 
tance of an STRE from the edge of a nucleosome is consis- 
tent with partial unwrapping of DNA from a nucleosome 
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in vivo. Finally, we note that Msn2 fails to bind to STREs 
lying in NDRs of limited width. This may suggest that such 
NDRs may be stably obstructed by factors other than nu- 
cleosomes or by nucleosomal subparticles of non-canonical 
length. In short, well-positioned nucleosomes, as well as 
other unidentified chromatin features, occlude STRE bind- 
ing of Msn2 and thereby dictate which of the many elements 
are available for promoting regulatory changes in response 
to stress-induced Msn2 entry into the nucleus. 

While nucleosomes generally restrict access of Msn2 to its 
cognate sites on the genome, Msn2 can serve to alter nucle- 
osome positioning. Shivaswamy et al. (11) showed that nu- 
cleosomes over some STREs were lost following heat shock 
and suggested that such loss allowed binding of Msn2 and 
subsequent transcriptional activation, arguing for a passive 
role for Msn2 in nucleosome remodeling. However, by ex- 
amining nucleosome positioning in both a wild-type and an 
msn2 msn4 strain, we have shown that in many cases nu- 
cleosome repositioning is dependent on Msn function and 
thus that Msn2, Msn4 or both play an active role in nucle- 
osome remodeling. Our results are consistent with those of 
Huebert et al. (41) who found that Msn2 binding to STREs 
occluded by nucleosome preceded nucleosome removal fol- 
lowing peroxide treatment. We find that Msn2 does not bind 
to STREs located under the central core of nucleosomes but 
rather can bind to sites located near the periphery of the 
nucleosomes. This suggests a model in which partial un- 
wrapping of nucleosomes in vivo allows initial binding of 
Msn2 to STREs at the edges of nucleosomes. Subsequent 
recruitment of chromatin modifying and remodeling pro- 
teins by Msn2 would lead to reposition or eviction of the oc- 
cluding nucleosome. In a similar way, recruitment of Msn2 
to STREs in NDRs could result in restructuring the local 
chromatin environment. Thus, while nucleosomes obstruct 
Msn2 binding, Msn2 binding can evict nucleosomes from 
some STREs and can alter the nucleosome positioning in 
its vicinity. 



Global measurement versus individual responses 

The studies described in this report measure Msn2 binding 
and nucleosome positioning over the entire population of 
cells in a culture and thus provide a measure only of the av- 
erage behavior over all cells. However, microfluidics-based 
studies demonstrate that stress elicits a variable Msn2 re- 
sponse, with some cells exhibiting substantial and sustained 
Msn2 nuclear occupancy while others showing no accu- 
mulation (34,35). Moreover, Msn2 shows bursting behavior 
in which the cellular cohort of Msn2 makes random tran- 
sient forays into the nucleus. Thus, the pattern of binding 
we observe by ChlP-Seq likely does not reflect the pattern 
of Msn2 genome occupancy in any individual cell, partic- 
ularly since the number of Msn2 molecules is apparently 
less than the number of bound sites identified in our and 
other studies. This raises the possibility that the cohort of 
genes induced byMsn2 may differ significantly among dif- 
ferent cells, either in a completely random fashion or in dis- 
tinct subsets of responses. Such diversity in transcriptional 
output may allow genetically identical cells in a popula- 
tion to respond individually and in quite distinct ways to a 



common stress, enhancing the number of survival strategies 
available to the population of cells as a whole. 
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